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ABSTRACT 

(j, Columbae is a prototypical weak-wind O-star for which we have obtained a high-resolution X-ray spectrum 
with the Chandra LETG/ACIS instrument and a low resolution spectrum with Suzaku. This allows us, for 
the first time, to investigate the role of X-rays on the wind structure in a bona fide weak-wind system and 
to determine whether there actually is a massive, hot wind. The X-ray emission measure indicates that the 
outflow is an order of magnitude greater than that derived from U V lines and is commensurate with the nominal 
wind-luminosity relationship for O-stars. Therefore, the "weak-wind problem" — identified from cool wind 
UV/optical spectra — is largely resolved by accounting for the hot wind seen in X-rays. From X-ray line profiles, 
Doppler shifts, and relative strengths, we find that this weak-wind star is typical of other late O dwarfs. The 
X-ray spectra do not suggest a magnetically confined plasma — the spectrum is soft and lines are broadened; 
Suzaku spectra confirm the lack of emission above 2 keV. Nor do the relative line shifts and widths suggest 
any wind decoupling by ions. The He-like triplets indicate that the bulk of the X-ray emission is formed rather 
close to the star, within 5 stellar radii. Our results challenge the idea that some OB stars are "weak-wind" stars 
that deviate from the standard wind-luminosity relationship. The wind is not weak, but it is hot and its bulk is 
only detectable in X-rays. 

Subject headings: stars: mass-loss — stars: early-type — stars: individual (p Col) — X-rays: stars 



1. INTRODUCTION 

The outflow of stellar winds from massive OB-type stars 
is an important process which affects both the chemical 
enrichment and kinetics of the interstellar medium (e.g. 
iLeitherer. Robert & Drissen| [r992). The mass-loss itself is 
enough to change the evolution of the star, which ends its 
life in a supernova explosion, also profoundly changing its 
environment. Hence, quantitative understanding massive star 
winds is important not only as a basic component of stellar 
astrophysics, but also for understanding cosmic feedback on 
galactic scales throughout cosmic history. 

While some basic physics of stellar winds in massive stars 
is well established — that winds are accelerated by photoelec- 
tric absorption of the intense ultraviolet radiation field by a 
multitude of metal lines, that an instability can lead to wind- 
shocks which generate X-rays — there are still puzzles to be 
solved. One of these is the "weak-wind" problem, in which 
UV line diagnostics clearly show a wind signature in clas- 
sical P Cygni line profiles, but modeled mass loss rates can 
be discrepant by more than an order of magnitude from val- 
ues expected based on O-star statistical trends and theoretical 
foundation s, specifically the wind momentum-luminosity re- 
lation (e.g. iPuls et all 1 19961) . Factors of a few in mass loss 
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are enough to be significant for stellar ev olution and cosmic 
feedback (e.g. lPuls. Vink & Naiarroll2008l) . 

It has long been known that photoionization by X-rays 
can alter the ionization b alance in the w i nd regions where 
the U V lines are formed (IWaldronl 119841: iMacFarlane et alJ 
fl99"l . There is a theoretical degeneracy in that different val- 
ues of mass loss rate (M) and X-ray luminosity (L x ) can 
produce very similar UV line profiles (IPuls. Vink & Najarrol 
I2008t iMarcolino etaD 120091) : direct knowledge of the X- 
ray spectrum is thus important for reliable determination 
of wind parameters. Another possibility is that cool and 
hot plasma emission originate from different volumes or 
densities; that is, clumping can affe ct the interpretation 
(lHamann. Feldmeier & Oskinoval2008[). ^ C ol belongs to the 
weak- wind domain defined by iLucyl (120101) in which a star's 
rate of mechanical energy loss in a radiatively driven wind 
is less than the radiative output from nuclear burning; Lucy 
showed that there is a huge disparity between the theoretically 
exp ected M and values derived from UV and optical spec- 
tra. iLucvl ( 120121) developed a phenomenological model and 
suggested that in low-luminosity O-type stars, the volumetric 
roles of hot and cool gas are possibly reversed compared to 
O-type stars of high luminosity; thus in the weak-wind stars, 
a larger volume is occupied by the hot gas than by the cool 
gas. 

/i Col (HD 38666) is an 09.5 V runaway (and single) star 
and is one of the we akest wind Galactic O-type stars (e.g., 
see Figs. 39 and 41 in lMartins et al.ll2~005l) . These factors are 
what motivated our spectroscopic study of the prototypical 
weak-wind system, /xCol, at high-resolution with Chandra 
and at low resolution but greater sensitivity at higher energies 
with Suzaku. In this letter, we concentrate on the primary 
empirical results from the X-ray spectral analysis of ^iCol. A 
subsequent paper will investigate the influence of the X-rays 
on the cool wind component (Todt et. al, in preparation). 
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Figure 1. /iCol C/icmrfra/LETG/ACIS spectrum (left); black: observed count-rate; thin-gray line (red in the online version): model; the lower panel shows 
the x 2 residuals against a broken powerlaw APEC model modified for photoexcitation of triplets. Lines are broader than the instrumental width. The Suzaku 
spectrum is shown on the right (black). For conciseness, we have summed the counts from the three detectors; such is not recommended when fitting, but it 
provides a good summary visualization of the data. The folded C/ianc/ra-derived model is in gray (red in the online version), and residuals below. We emphasize 
that the C/icffirfra-derived model was not fit to the Suzaku spectrum, only folded through the response to provide model counts and residuals; the model agrees 
very well without any adjustments. It is significant that there is little or no flux detected above 2 keV where Suzaku has substantial sensitivity. 



Table 1 

H Col Properties 



Property 



Value 



Property 



Value 



Spectral Type 


09.5 V 


d [pc] a 


408 


T eff [K] b 


33000 


R/Rq^ 


4.6 


Iog(L bol /i )<v 


4.4 


tioo[kms 1 ] b 


1200 


M uv [M /yr] b 


10 -9.5 


JV H [cm~ 2 ] c 


5 x 10 19 


/4l-40A)[cgs] d - f 


6 x 10 -13 


log(L x /L bol ) d 


-6.9 


10 54.1 


T max [MK] d 


4.4 




1600 (±275) 






f/i{0 vn)8 


< 0.01 


r(0 VIl)8 


< 3.3 


//zfNe ix)s 


0.04-0.14 


r(Ne ix)8 
r(Mg Xl) g 


2.3-4.4 


//i(Mg xi)« 


> 0.2 


> 2.1 



a From the Hipparcos parallax, as re-evaluated by van Leeuwen] <2007l) 

b |Maninsetalji200l) 

c ICassinelli et an Tl 994); How k, Savage & Fabianl (T553) 
d This work 

c Adjusted for the adopted distance. 

f The flux is as observed at Earth, with foreground absorption. 
g f/i gives the ratio of the forbidden to intercombination line fluxes, and 
r is the radius of formation in units of the stellar radius, as derived from 
PoWR models. 

2. ANALYSIS 

We have obtained a 232 ks exposure of /iCol with the 
Chandra LETG/ACIS instrument (observation IDs 12349, 
12350, 13422, PI L. Oskinova). Figure Q] shows the count- 
rate spectrum. Since /i Col is a single star, there are no am- 
biguities present as when interpreting observations of binary 
systems with composite spectra or colliding wind emission. 
Relevant stellar properties are given in Table Q] 

There are several key X-ray spectral stellar wind diag- 
nostics. The line profile is sensitive to the wind opacity 
and v elocity field dMacFarlane et al.ll 199 lb lOwocki & Cohenl 
120011) : the line centroid and width are useful proxies, be- 
ing sensitive to wind parameters governing the detailed line 
shape. Emission line strengths are indicative of plasma tem- 
peratures and elemental abundances. The continuum at the 
shortest wavelengths available (2-10 A, 1-6 keV) is also very 
sensitive to the highest temperatures present. The He-like 
triplets are sensitive to electron density and the UV radiation 
field through collisional and photoelectric excitation which 
can depopulate the forbidden-line level, w eakening it while 
strengthening the intercombination lines dGabriel & Jordan! 
1 19691 iBlumenthal. Drake & Tucker! 1 19721) : in O-stars, the 
UV field typically dominates the depopulation and hence 



He-like line ratios are dia gnostic of radius of formation 
(IWaldron & Cassinellill200ll) . Very close to the photosphere, 
density effects could also become significant. 

We have fit a global model to the X-ra y spectrum, us- 
ing s tandard products produced by CIAO dFruscione et all 
120061 ver sion 4.3, and associated calibration database), us- 
ing ISIS dHouck & Denicoial 120001) and the collisional ion- 
ization equilibrium emiss ivities in AtomDB (version 2.0, 
ISmith & Brickhousell2008l) . For the plasma model, we used 
a broken powerlaw emission measure distribution (EMD; 
which in differential form is defined as n e rihdV / dT), with 
variable abundances for significant ions, a Doppler shift, and 
line profiles defined by a global Gaussian Doppler broaden- 
ing term. Our models show that the absorption of X-rays in 
the cool stellar wind is negligible, and spectra can be well 
fitted neglecting wind absorption. Abundances, velocity, and 
broadening were common over all temperature components. 
The resulting fit is shown as the red curve in Figure Q] The 
emission measure distribution is shown in Figure [2] 

Abundances (refer enced to solar photospheric values of 
lAsplund etaf] (120091) ) were about half of solar (O, Si, and 
Fe) or near solar (N, Ne) — we could not obtain a good fit 
using abundances all set to solar. This may be due to the 
adopted smooth functional form of the emission measure dis- 
tribution, since the EMD and abundances are somewhat de- 
generate. Trial fits with discrete temperature components 
were also poorer when using solar abundances. To explore 
this somewhat further, since the integrated emission measure 
is fundamental to our main result, we evaluated the 90% con- 
fidence limits of the model normalization, froze this at each 
of the high and low limits, and re-fit the spectrum to ob- 
tain new EMD and abundances. We could have half the 
integrated EM with relative abundances increased to 0.7- 
1.5, or we could have about triple the best-fit EM, also with 
slig htly modified abundances. ( We note that from UV spec- 
tra, iFitzpatrick & Mas"sal d 19991) also found low abundances 
for /i Col, but did not believe their results plausible and at- 
tributed them to model deficiencies.) 

The broken powerlaw model is empirically justified in that 
it provides a necessary multi-thermal model and does fit 
well with relatively few parameters. Such an EMD can be 
physically justified by hydrodynamic shock models which 
predict a wide range of temperatures over a large range 
of radii dFeldmeier. Puis & Pauldrachl 119971) . Other OB- 
stars have also shown similar, empirically determined EMD 
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Figure 2. The best fit emission measure distribution is given by the thick 
solid line. Circles mark the 1 dex temperature intervals, the resolution of 
the emmissivity database, and the value plotted is the emission measure inte- 
grated over 0. 1 dex. Vertical bars near the top mark the temperatures of max- 
imum emissivity for detected emission lines, and the abscissa's range spans 
the temperatures where these lines have greater then 50% of their maximum 
emissivity. The instruments have significant sensitivity from 3-6 A such that 
EM above 10 7 K would be apparent in lines and continuum. 



dWoidowski & Schulzll2005l) . The important point here is to 
obtain an order-of-magnitude estimate of the EM, and that can 
be done with a variety of plausible models. Also given that 
the X-ray emitting plasma is optically thin, it does not mat- 
ter (for the EMD) where the emission originates, under the 
assumption that abundances are independent of temperature 
and density. 

Model properties are given in Table [T] the volume emission 
measure, EM X , is integrated over the temperature range from 
1-100 MK. 

We have observed /.t Col using the Suzaku satellite for 26 ks 
(observation ID 405059010, PI L. Oskinova). We used the 
HEASoft 6.9 and the CALDB xis20100123 for the XIS anal- 
ysis; the HXD did not detect any signal. The Suzaku spectrum 
is consistent with the above model, with no flux above 2 kcV. 
Figure [TJ shows the Suzaku spectrum, the folded Chandra- 
derived model, and residuals; there is good agreement. 

Global fits cannot necessarily provide detailed information 
present in individual lines, particularly if lines from different 
ions have different characteristics from local conditions (e.g., 
velocity and temperature gradients). Hence, we also fit the 
stronger lines in the LETG/ACIS spectrum with Gaussians 
(folded through the instrument response), with a continuum 
derived from the global plasma model. The parameters rele- 
vant here are the centroids and widths (Figure[5]i. A non-zero 
centroid indicates a wind, being skewed to the blue by disk oc- 
cultation of the receding wind and by absorption in the wind. 

We als o computed a model line shape for an expand- 
ing wind (lOskinova. Feldmeier & Hamannll20 06) and fit this 
to the strongest, isolated feature in the spectrum, O VIII 
(18.967, 18.973 A) by adjusting the line position, flux, and 
scaling the width (Figure@]i. We used an un-clumped model 
profile with (3 = 0.7 and R /R± = 1.1. We find the wind 
is thin enough that a small asymmetry and shift is due to the 
stellar disk occultation of the receding hot wind. The blueshift 
in the Gaussian fit to O VIII of —84 km s _1 (Figure|3]l is con- 
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Figure 3. /i Col emission line Doppler offset and width; errorbars give 90% 
confidence limits. Due to lower resolution at shorter wavelengths, the Mg XI 
and Si XIII lines are poorly constrained. The systemic line-of-sight velocity 
of 109.2 km s _1 (Evans 1967) has been removed. 

sistent with the wind profile model, whose best fit is at the 
expected line position (i.e., no offset, within one standard de- 
viation accuracy of 50 km s _1 ). We did not achieve a good fit 
for the UV-derived v x = 1200 km s -1 , but had to increase 
the hot wind velocity to 1600 (±275) kms^ 1 . UV spectra do 
not allow precise determin ation of ; publishe d values range 
from 1000-2000 km s" 1 dMartins et al.l 120051) . Our PoWR 
models showed that the C IV, Si IV, and N V lines can be 
equally well fit with Vqo of 1200 km s _1 or 1600 km s _1 . We 
can also achieve an equally good fit to O VIII for j3 = 1, but 
with = 2800 km s -1 . More detailed analysis with more 
lines and a variety of model assumptions are required. For a 
conservative approach, we prefer the lower v^. 
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Figure 4. /i Col O VIII profile (black), wind-profile fit (gray; red in the online 
version) and residuals (lower panel). The model has a line center offset of 
24 (±54) km s" 1 , and Voo = 1600 (±275) km s _1 . The model profile 
was for = 0.7, Ro/R* = 1.1, with a smooth, un-clumped wind. 

The final line measurements of use are the He-like triplet 
ratios. We have measured these for O VII, Ne IX, and Mg XI 
(Si XIII is too weakly exposed and significantly blended to 
be useful). The forbidden lines are much weakened: The 
O VII forbidden line was not detected, and the Ne IX forbid- 
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den line was very weak. In both cases, the intercombination 
lines are as strong as the resonance lines. The forbidden-to- 
intercombination ratios (f/i) provide limits on the radii of 
formation. We applied a simple model assuming triplet for- 
mation occurs at a single radius in the wind, using only ex- 
citation by photospheric UV flu xes from model atmospheres 
such as TLUS TY or ATLAS-9 dLanz & Hubeny||2007l 120031: 
lKuruczlfT979l) . We als o used more sophisticated treatment 
using the PoWR code (lHamann & Graf ener 2004) in which 
the wind model calculates ionization structure, solves the en- 
ergy equation, includes X-rays, includes photospheric emis- 
sion corrected for absorption lines, line blanketing, non-LTE 
populations, stellar disk limb darkening, and the diffuse UV 
radiation field. Each of the methods predicts the / / i -ratio vs 
radius r of formation. Values (which were nearly identical for 
all methods) are given in Table[T]with r, in units of the stellar 
radius. 

There is some temperature dependence in the triplet ratios, 
but it is below 10%, less than our measurement uncertainty. 
The point-formation model also ignores the effect of any dis- 
tributed emission. If there is contribution from further out in 
the wind, our / / i-derived radii will be overestimates (since / 
increases and i decreases at larger radii — the distributed emis- 
sion ratio would be larger than at any point at smaller radius). 

We have also looked for variability. The count rate in dis- 
persed photons appears constant; the rate cannnot change by 
more than 5-10% in time intervals of 1 ks or more. 

3. "WEAK- WIND" STARS ARE NOT LOW MASS-LOSS WINDS 

The He-like triplets place the X-ray emission within about 
2-5 stellar radii of the photosphere. Line widths are con- 
sistent with X-rays being formed at small radii. We as- 
sume that the hot plasma obeys the continuity equation, that 
the standard velocity law holds for the hot plasma, namely 
v x = Voc,x(l — bR^/r) 13 , and that the hot plasma exists only 
above some inner radius Rq > For Voo = 1200kms _1 , 
b = 0.97 (a value based on the ratio of the sound speed at i?* 
to Voo, which provides a non-zero wind velocity, but which 
is somewhat uncertain.) and j3 = 1, we expect wind veloc- 
ities in this region to be about 800-1300 km s -1 . Figure [3] 
shows the widths to be somewhat smaller at 300-700 km s _1 . 
This suggests that much of the X-ray emission originates very 
close to the star within the radiation-driven wind acceleration 
zone, or that the hot plasma does not follow the typical ve- 
locity law. The hot plasma, however, does expand with rela- 
tively high velocity since the X-ray emission line profiles are 
resolved. Detailed line-profile fitting using wind-models will 
be required to determine the structure in more detail. 

With the assumptions above and the definition of the X-ray 
emission measure (and /3 = 1 to allow an analytic integra- 
tion), we can derive a simple expression for the X-ray inferred 
mass loss rate in solar masses per year, also assuming that the 
hot wind is unclumped, and that the cool wind is optically thin 
to X-rays (as justified by the line profiles): 



M x = 3 x 10 9 Voo , x 



Rq \R+ 



1/2 



(1) 



in which Voo^ x is the hot wind's terminal velocity as de- 
termined from X-ray emission line profiles (in units of 
1000 km s -1 ), R± is the stellar radius, EM X is the emission 
measure of the hot plasma (in units of 10 54 cm -3 ), and b is 
the unitless parameter from the wind velocity law. Using 
Rq = i?*, b = 0.97 and values from Table Q] we infer that 



M w 2 x 10~ 9 Mp, yr -1 . This is 6 times the value derived 
from the UV by Martins et al. ( 2005) (or more, accounting for 
subsequent revisions in distance and stellar radius). If the hot 
wind begins somewhere above the photosphere, then the in- 
ferred mass loss rate will be even larger (e.g., about 20 times 
fori?o/i?* = 1.5). 

The above M x uses our EM X ~ 10 54 . If were to assume 
the much lower M uv , we would infer a much smaller emis- 
sion measure for a spherically symmetric wind. This means 
that most of the wind is hot (or denser) than the UV-emitting 
plasma, a si tuation also noticed in the study of main sequence 
B stars (e.g. lCassinelli et al.|[l994l) . 

We conclude that the hot wind of [i Col must have a larger 
volume or greater density than the cool wind. The wind is not 
weak, but it is hot and its bulk is only detectable in X-rays. 

Our observations exclude other proposed weak-wind ex- 
planations, such as magnetically channeled wind shocks 
(MCWS) or frictionally decoupled winds. The X-ray spec- 
trum of /i Col is not characteristic of MCWS in which hot 
plasma is held in stationary structures close to the stellar 
photosphere and which reach high te mperatures from col- 
lision of funneled high velocity winds (|Babel &Montmerlel 
U997tlTownsend. Owocki & Ud-D oula 2007|). The pi asma of 
^iCol is not extremely hot (T > 10 MK) nor are lines un- 
shifted and unresolved. Decoupling of ions from neutrals 
could occur at low densities or very low metallicities, cre- 
ating a two-fluid system (c.f..lSpringmann & Pauldrachl 19921: 
lKrti5ka & KubatU200it iMartins et al.ll2004l) . Frictional heat- 
ing by decoupled ions seems unlikely: the plasma is not ex- 
tremely cool (< 1 MK) with different widths and velocities 
for different ions (within our sensitivity; Fig. |3). Given the 
X-ray EM X and likely radii of formation, the density is not as 
low as once presumed. The relative Fe abundance from our 
fi ts is about 0.5 solar, not low enough to cause decoupling. 

iLucvl (1201 Oh reduced the weak-wind discrepancy through 
theoretical arguments. The X-ray spectra of /xCol support 
this with empirical evidence for a dominant hot wind, even 
given some uncertainty from still poorly determined factors of 
clumping, wind veloci t y law, and location of the wind base. 
iDrew. Hoare & Denbyl (119941) suggested that X-rays may in- 
crease the wind ionization at the critical point, lowering the 
effective radiative acceleration and mass-loss rate. Our Chan- 
dra observations resolve the broad X-ray emission lines, indi- 
cating that the hot plasma expands at high velocities, perhaps 
even exceeding the cool wind velocity determined from UV 
spectra. 

Ou r empi rical findings, reached in ignorance of work by 
ILucvl d2012h . are in line with his claim that in late-type O- 
dwarfs the ambient wind is heated to temperatures of few MK 
at radii > 1.4 i?*, with cool radiatively-driven gas being con- 
fined to dense clumps and small volume filling factor. Fur- 
ther out in the wind in his model, cool clumps are destroyed 
by heat conduction from the hot plasma, and the outflow is 
dominated by a hot thermal wind which reaches a supersonic 
terminal velocity of ~ lOOOkms -1 . 

In future work, we will refine our order-of-magnitude esti- 
mates using detailed wind models. For a relatively thin wind, 
we do not expect significant qualitative changes. Some ob- 
servational uncertainties need to be refined: we only have 
both upper and lower limits for the formation radii of Ne IX; 
whether other species form at 5 i?* or below 2 i?* is important 
for wind shock models. Higher spectral resolution would be 
of use, especially for the hottest He-like triplet in the spec- 
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tram, Si XIII, to probe the deepest layers. Higher resolution 
would better determine line shifts and widths, and better con- 
st rain wind structure. 

iNajarro. Hanson & Pulsl (1201 ll) recently classified crOri 
AB (HD 37468, 09.5 V + B0.5 V) as a weak-wind sys- 
tem, based on infrared spectra, revising its mass-loss rate 
do wnward by more t h an tw o orders of magnitude from that 
in iHowarth & Prinjal dl989l) . This st ar has also been ob- 
served at high resolution in X-rays dWaldron & Cassinellil 
I2007t IZhekov & Pallall2007b ISkinner et al.ll2008l) . Based on 
these prior works and our own analysis of the Chandra/HETG 
spectrum, we find that a Ori has X-ray spectral characteristics 
very similar to those of /i Col. This further corroborates our 
conclusions that the weak-wind phenomenon is due to a prop- 
erty of the cool plasma being a minor constituent of the wind. 

4. CONCLUSIONS 

Our results challenge the idea that some OB stars are 
"weak-wind" stars that deviate from the standard wind- 
luminosity relationship. From high-resolution X-ray spec- 
trum of yu Col, specifically He-like lines and the total emission 
measure, this star does not appear unusual relative to other O- 
stars except for its weak-wind status. Its X-ray emission mea- 
sure, line widths, and centroids are in good agreement with 
the OB main sequence star results of I Waldron & Cassinellil 
(120071) . The volume emission measure of the X-ray emitting 
plasma must be very much larger than the cool, UV-emitting 
plasma, and we believe that the weak-wind problem is re- 
duced or eliminated when the hot and dominant component 
of the wind is taken into account. The wind is not weak, but 
it is hot and its bulk is only detectable in X-rays. 
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